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Abstract One of the most dramatic events in the life 
of a low-mass star is the He flash, which takes place at 
the tip of the red giant branch (RGB) and is followed by 
a series of secondary flashes before the star settles into 
the zero-age horizontal branch (ZAHB). Yet, no stars 
have been positively identified in this key evolutionary 
phase, mainly for two reasons: first, this pre-ZAHB 
phase is very short compared to other major evolution- 
ary phases in the life of a star; and second, these pre- 
ZAHB stars are expected to overlap the loci occupied by 
asymptotic giant branch (AGB), HB and RGB stars ob- 
served in the color-magnitude diagram (CMD). We in- 
vestigate the possibility of detecting these stars through 
stellar pulsations, since some of them are expected to 
rapidly cross the Ccpheid/RR Lyrae instability strip in 
their route from the RGB tip to the ZAHB, thus becom- 
ing pulsating stars along the way. As a consequence of 
their very high evolutionary speed, some of these stars 
may present anomalously large period change rates. We 
constructed an extensive grid of stellar models and pro- 
duced pre-ZAHB Monte Carlo simulations appropriate 
for the case of the Galactic globular cluster M3 (NGC 
5272), where a number of RR Lyrae stars with high pe- 
riod change rates are found. Our results suggest that 
some - but certainly not all - of the RR Lyrae stars in 
M3 with large period change rates are in fact pre-ZAHB 
pulsators. 
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1 Introduction 



stars: Horizontal 



Low-mass stars, after exhausting their central hydrogen 
supply, develop increasingly more massive helium cores 
due to the continued supply of He-rich material from 
the H-burning shell that surrounds it. Eventually, the 
conditions in the core become ripe for the ignition of 
the He-burning triple-alpha reactions. Since the core 
at this point is electron-degenerate, the process does 
not take place quiescently, but rather in the form of a 
thermonuclear runaway, called He flash. In this process, 
the star will typically change its luminosity by several 
orders of magnitude; the surface temperature may also 
increase very significantly. The inner structure of the 
star undergoes dramatic changes in the process as well. 
And yet, in spite of the importance of this evolution- 
ary phase, to the best of our knowledge not a single 
such star has positively been identified to date. The 
main goal of the present study is, accordingly, to per- 
form the first systematic study of the expected prop- 
erties of pre-ZAHB stars, with a view towards their 
detection using a combination of pulsation and color- 
magnitu de diagram (CMD) prope rties. For more de- 
tails, see lSilva Aguirre et alj (|2008l ). 



2 Stellar Grid and Population Synthesis Code 



We u sed the Garching Stellar Evolution Code ( Weiss fc Schlattlfr 
20081) to construct the evolutionary tracks required to 
populate the pre-ZAHB phase. We adopted the mix- 
ing length theory with a value «mlt = 1-75, and a 
chemical composition of Z = 0.001 and Y = 0.25, as 
appropriate for M3. Our models were computed with 



2 



an initial mass on the zero-age main sequence (ZAMS) 
of 0.85 M Q and application of mass loss in its most 
widel y used formulation, namely the Reimers ( 19751 : 
19771) formula, 



<9M n L 



(i) 



where L is the stellar luminosity (in solar units), g the 
surface gravity (in cgs units), and R the surface radius. 
A thin interpolation grid was created by evolving stars 
from the ZAMS to the ZAHB for a range of ?7r, from 
0.0 to 0.60. We adopted a Gaussi an-shaped mass dis- 
tribu tion at the tip of the RGB (jValcarce fc Catelan 
20081) ). with mean value < M > = 0.647 M and a = 
0.022 M Q . We generated the masses for the stars pop- 
ulating the prc-ZAHB randomly, and interpolated over 
the tracks to obtain their luminosity and temperature. 
We then computed the blue ed ge of the instability strip 



Caputo et al.l (|l987T ). after applying a correc- 
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tion of -200 K. The width of the IS has been taken as 
Alog(T c ff) = 0.075. Periods and period change rates 
are calculated for the stars contained withi n the IS us- 
ing the formulation of ICaputo et al.l (|1998l ). In Fig. [1] 
we present the result of one of our simulations with one 
of our evolutionary tracks overplotted. 




Fig. 1 Results for one simulation with 457 pre-ZAHB 
stars. Red triangles indicate our simulated stars, the solid 
line depicts the evolutionary track for rjn = 0.41, and the 
green dotted line shows the position of the ZAHB for differ- 
ent mass values. Blue dashed lines are the instability strip 
boundaries (schematic). 



3 Empirical Data 

Magnitudes and c olors for the variabl e stars in M3 
were obtained from lCacciari et aL ( 2005 ). For the non- 
variable stars, we use d the photometry provid ed by 
Buonanno et all (jl994h and iFerraro et al.l (jl997h . M3 
has ~530 HB stars whose position in the CMD have 
been measured. If we consider lifetimes, pre-ZAHB evo- 
lution lasts ~1.4 Myr, while the HB phase takes ~80 
Myr. Thus, one should expect to find one pre-ZAHB 
star for every 60 or so bona fide HB stars. 



4 Results 

In Fig. [5] we present a histogram that shows the ex- 
pected period change rate distribution for an M3-like 
input mass distribution. As can be seen, most of the 
pre-ZAHB stars have period change rates contained in 
the [-1,0] d/Myr range, meaning that evolution occurs 
mainly towards higher temperatures and/or decreasing 
luminosity, as expected for these types of stars. 
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Fig. 2 Distribution of the period change rates. The middle 
panel zooms in around the peak of the distribution, while 
the lower panel shows the wing of the distribution for the 
more extreme (negative) values. 

Therefore, among 530 HB stars in M3, approxi- 
mately 9 should be in the pre-ZAHB phase - two of 
which being pulsating (RR Lyrae-like) stars. We em- 
phasize the fact that AP/At absolute values higher 
than approximately 0.1 d/Myr arc completely unex- 
pected on the basis of canonical evolution, but have 
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often been observed among field and cluster RR Lyrae 
stars. 

In Fig. [31 the time distribution along the evolution- 
ary track corresponding closely to the most likely mass 
among M3's HB stars is presented. A red dot was 
placed at every 0.01 Myr of evolution, starting from 
the He flash down to the ZAHB. 





Fig. 4 Distribution of the M3 RR Lyrae stars on the CMD. 
Filled red circles: stars with AP/At < —0.15 d/Myr; open 
blue circles: stars with AP/At > +0.15 d/Myr; triangles: 
RR Lyrae stars with smaller (or undetermined) AP/At val- 



5 Conclusions 



Fig. 3 Pre-ZAHB evolution on the Hertzsprung-Russell 
diagram for the 77R. = 0.41 evolutionary track. Red dots are 
plotted for every 0.01 Myr of evolution since the main flash. 

As can be seen, the most extreme negative values 
of the period change rate can be obtained during the 
(very) fast evolution along the secondary loops. How- 
ever, this is a relatively unlikely situation, since pre- 
ZAHB stars spend most of their time in the blueward 
evolution that takes place after the last secondary flash 
(just prior to the ZAHB proper). This implies that 
most of the pre-ZAHB variables will have negative pe- 
riod change rates in the aforementioned [-1,0] d/Myr 
range. 

In Fig. IU we present the RR Lyrae stars of M3 
with three of our evolutionary tracks ovcrplotted, whose 
associated mass values at the tip of the RGB make them 
representative of the majority of the pre-ZAHB stars 
in the cluster. As can be seen, the positions of these 
variables in the CMD, as well as their period change 
rate values, are broadly consistent with at least some 
of them being pre-ZAHB stars. Indeed, according to 
our simulations, at least ^2 of them very likely are pre- 
ZAHB pulsators. 



Our pre-ZAHB Monte Carlo simulations show that, 
for an M3-like HB mass distribution, the pre-ZAHB 
CMD positions should be basically indistinguishable 
from those of bona fide HB stars in the cluster. We 
predict that, in the specific case of M3, ~9 pre-ZAHB 
stars should be present, ~2 of which most likely are 
disguised as RR Lyrae stars with relatively large (and 
negative) period change rates. 

Acknowledgements Support for M.C. is provided 
by Proyccto Basal PFB-06/2007, by FONDAP Ccntro 
de Astrofisica 15010003, and by Proyecto FONDECYT 
Regular #1071002. 



•1 



References 

Buonanno, R., Corsi, C.E., Buzzoni, A. et al. 1994, Astron. 

Astrophys., 290, 69 
Cacciari, C, Corwin, T.M., & Carney, B.W. 2005, As- 
tron. J., 129, 267 
Caputo, F., de Stefanis, P., Paez, E.,& Quarta, M.L. 1987, 

Astron. Astrophys., 68, 119 
Caputo, F., Marconi, M., Santolamazza, P., 1998, Mon. Not. 

R. Astron. Soc, 293, 364 
Ferraro, F.R., Carretta, E., Corsi, C.E. et al. 1997, Astron. 

Astrophys., 320, 757 
Reimers, D. 1975, Mem. Soc. R. Liege 6. Ser., 8, 369 
Reimers, D. 1977, Astron. Astrophys., 57, 395 
Silva Aguirre, V., Catelan, M., Weiss, & A., Valcarce, A.A.R 

2008, Astron. Astrophys., 489, 1201 
Valcarce, A.A.R, Catelan, M. 2008, Astron. Astrophys., 487, 

185 

Weiss, A., Schlattl, H. 2008, Astrophys. Space Sci., 316, 99 



This manuscript was prepared with the AAS IATf^X macros v5.2. 



